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^ . Abstract. A number of problems of solar and stellar dynamo theory 

O \ are briefly reviewed and the current status of possible solutions is dis- 

^ | . cussed. Results of direct numerical simulations are described in view 

of mean-field dynamo theory and the relation between the a-effect and 
■ the inverse cascade of magnetic helicity is highlighted. The possibility 

of 'catastrophic' quenching of the a-effect is explained in terms of the 
constraint placed by the conservation of magnetic helicity. 
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1. Introduction 

O . 

The solar dynamo problem has always been plagued by a number of problems. 
One of the most outstanding is the question of why the solar dynamo wave trav- 
els toward the equator and not the other way around. According to ail-dynamo 
theory the dynamo wave travels equatorward only if the sign of the product of 
alpha-effect and radial angular velocity gradient, d£l/dr, is negative. However, 
cyclonic convection gives rise to a positive a-effect (Parker 1955, Steenbeck, 
Krause, & Radler 1966), so one needs a negative sign of dfl/dr. Thus, a prob- 
lem has arisen since the mid-eighties (Parker 1987), i.e. since helioseismology 
took away the freedom of dynamo theorists to adopt a suitable rotation profile. 
This left only the choice of adopting an appropriate profile (and sign!) for the 
alpha-effect. Not surprisingly, dynamo theorists were labelled as being able to 
reproduce almost anything, but in reality they were not even able to do that, 
because there was one more problem: the phase relation between poloidal and 
toroidal fields (Stix 1976). The sign of the mean radial magnetic field, B r , is 
observed to be almost always opposite to the sign of the mean toroidal field, 
Ba,- The phase relation between B r and B^ is determined directly by dCl/dr, 
because the shear turns radial field into toroidal. This problem is actually rather 
general in its nature and quite independent of dynamo theory. 

One possible solution to this dynamo dilemma (Parker 1987) may lie in the 
possibility that the sense of the dynamo wave can be determined by the sense 
of the meridional circulation: a poleward circulation in the upper layers corre- 
sponds to equatorward motion in the lower layers, and if the sunspot activity 
is governed by fields at the bottom of the convection zone this circulation may 
control the migration of activity regions Durney 1996, Choudhuri et al. 1996). 
A systematic survey of solutions with meridional circulation (Kiiker et al. 2000) 
has revealed that the sense of the dynamo wave is reversed only if the following 
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conditions are met: the magnetic eddy diffusivity is small enough, the circula- 
tion is strong enough, and the a-effect is sufficiently supercritical. If this is the 
case, Kiiker et al. (2000) also find B r and to be approximately in antiphase, 
as observed. 

Over the past ten years some rather more fundamental problems have shad- 
owed the attempts to model closely the solar field geometry: does the a-effect 
really work in high Reynolds number flows like the sun, or is it totally suppressed 
by the magnetic field? 



2. The helicity constraint 

A serious problem that has only recently received attention results from magnetic 
helicity conservation. The problem is relatively easily explained and applies 
to all fields that have magnetic helicity. In particular, it applies to the large 
scale magnetic fields generated by a 2 or afi dynamos, but not, for example, to 
fields generated by a small scale dynamo which may produce significant helicity 
fluctuations, but no net magnetic helicity. 

Magnetic helicity can only change resistively or through flux on the bound- 
aries. Thus, a large scale field with magnetic helicity can only be generated 
on a resistive timescale or, if there are suitable losses of the right kind on the 
boundaries, field generation can proceed on a dynamical timescale. The latter 
alternative is not straightforward, because such losses always involve losses of 
magnetic energy too. This latter alternative will be discussed later. First, how- 
ever, we shall look more quantitatively at the consequences of resistively slow 
changes of helicity in closed or periodic boxes. 

The condition of helicity conservation together with the assumption that 
the large scale field is helical result in a simple condition for the energy of the 
mean magnetic field (Brandenburg 2000, hereafter referred to as B2000) 

(B 2 )/B 2 cq « k f e AB [l - exp(-27 ? A ; 2 (t - t sat )] , (1) 

where k{ is the wavenumber of the (kinetic) energy carrying eddies (e.g. the forc- 
ing wavenumber), k\ is the smallest possible wavenumber, rj is the microscopic 
magnetic diffusivity, t is time, t sa t is the time when the field at small and inter- 
mediate scales saturates, and Iab quantifies the degree to with the large scale 
field is helical. This relation is rather general and independent of the actual 
model of field amplification. If the field is not fully helical, then Iab will be 
reduced, but the important point here is that full saturation is only obtained 
after a large scale ohmic diffusion time, l/r/kf. 

The important point is that all box simulations with helicity have to obey 
this behavior. This is shown in Fig. |l| for one particular example, but this 
constraint is indeed very well obeyed by all other simulations presented in B2000. 

The helicity constraint has two properties: full saturation of the large scale 
field (as opposed to the small scale field) happens on a resistive time scale, 
I/77 A; 2 , and secondly, the energy of the maximum field strength in units of the 
equipartition field strength is given by the product kflAB, and for maximally 
helical fields we have £ab = k^ 1 . In our particular runs we have ki/k\ = 5, so we 
can reach maximally five times supercritical energies of the large scale field. (In 
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Figure 1. Evolution of (B ) for Run 3 of B2000, compared with 
the constraint imposed by helicity conservation (dashed lines). Also 
plotted is the result of a mean-field dynamo with 'catastrophic' alpha 
and eta quenching (dotted lines). 



B2000 there is also run with kf/k% = 30, which does indeed give correspondingly 
higher field strengths.) The reason we have this factor 5 here can be motivated 
as follows. The large scale dynamo is driven by some effective a which has 
contributions from the kinetic helicity and the current helicity of the small scale 
field which combine to a residual helicity (Pouquet et al. 1976, Field et al. 
1999). Near equipartition the residual helicity is very small, so kinetic helicity 
and small scale current helicity nearly balance each other. The kinetic helicity, 
(u ■ u), is proportional to k{(u 2 ), so we know how big that is. On the other 
hand, the small scale current helicity, (j • b), is exactly the term which drives 
magnetic helicity, (A ■ B). In the steady state (A ■ B) can no longer change, so 
{j ■ b) must be balanced with the contribution from large scales, (J • B), which 

determine then the large scale field strength, (J ■ B) = k\(B ). This is - in 
words - the reason why superequipartition field strengths can be attained. (For 
mathematical details see B2000.) 

2.1. Catastrophic quenching 

Exactly the same resistive behavior as in Eq. ^ can be reproduced using the o? 
dynamo equation (e.g. Moffatt 1978) with simultaneous quenching of a and the 
turbulent magnetic diffusivity, r/t, in the form 

"o r] t0 

a = =2 ) Vt = =o > l 2 J 

l + a B B 2 /Bi q l + r, B B 2 /Bl q 

where we assume as = f]B- We focus here on the case of homogeneous turbulence 
in which case the large scale magnetic field is a force-free field whose energy 

2 

density is very nearly uniform, so B is just a function of time and we can 
obtain the solution B = B(t) in closed form (B2000). 

A natural by-product of this model is that the quenching coefficient as 
(which is assumed to be the same as t]b) must be proportional to the magnetic 
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Reynolds number! This type of quenching was suggested earlier on phenomeno- 
logical grounds by Cattaneo & Vainshtein (1991) and Vainshtein & Cattaneo 
(1992), but the conclusion was thought to be that for dynamically important 
field strengths (\B\ — > B cq ) the two turbulent transport coefficients would be 
so strongly quenched that such fields cannot be the result of a mean-field dy- 
namo. That was the reason that such quenching was referred to as 'catastrophic'. 
We now see however that there is nothing catastrophic about it and that field 
strength even in excess of B cq can be generated. 

The significance of these results is that they provide an excellent fit to 
the numerical simulations; see Fig. |l] where we present the evolution of (B ) for 
Run 3 of B2000. The dynamo equations with appropriate quenching expressions 
can therefore be used to extrapolate to astrophysical conditions. The time r eq 
required to convert the small scale field generated by the small scale dynamo 
to a large scale field increases linear with the magnetic Reynolds number, Rrn • 
Apart from some coefficients of order unity the ratio of r eq to the turnover time 
is therefore just R m . For the sun this ratio would be 10 8 — 10 10 . However, 
before interpreting this result further one really has to know whether or not 
the presence of open boundary conditions could alleviate the issue of very long 
timescales for the mean magnetic field. Furthermore, it is not clear whether 
the long timescales discussed above have any bearing on the cycle period in 
the case of oscillatory solutions. The reason this is not so clear is because for a 
cyclic dynamo the magnetic helicity in each hemisphere stays always of the same 
sign and is only slightly modulated. It is likely that this modulation pattern is 
advected precisely with the meridional circulation, in which case the helicity 
could be nearly perfectly conserved in a lagrangian frame. This would support 
the suggestion of Durney (1995) and Choudhuri, Schiissler, & Dikpati (1995) 
that the dynamo wave travels mainly because of meridional circulation. 

3. Quadratic or cubic quenching? 

When the first nonlinear mean-field dynamo models with a-quenching were cal- 
culated numerically (e.g., Jepps 1975, Ivanova & Ruzmaikin 1977), a simple 
quenching formula of the form a ~ 1/(1 + B 2 ) was used. This was really noth- 
ing more than just a simple fix to the problem that a quadratic formula of the 
form a ~ 1 — B 2 may overshoot and change sign. We emphasize that a quench- 
ing formula of the form a ~ 1/(1 + B 2 ) was never obtained rigorously. Instead, 
the correct quenching behavior, at least in the low magnetic Reynolds number 
limit, had cubic limiting behavior (Moffatt 1972, Riidiger 1974). 

We are now for the first time in a position to check whether the limiting 
behavior for higher magnetic Reynolds numbers is quadratic or cubic. In Fig. || 
we plot the evolution of the magnetic energy for a model with cubic a-quenching, 

a = a /[l + (a B B 2 /B 2 q ) n }, n = 3/2. (3) 

2 

It turns out that in this case the evolution of B does not agree with the lim- 
iting behavior enforced by the helicity constraint ([]]); see Fig. |2[ This confirms 
that the phenomenological expression (||) is actually correct, at least for the 
homogeneous a 2 dynamo. (We shall see later that for inhomogeneous dynamos 
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with shear the tensorial nature of a-effect and turbulent diffusivity cannot be 
ignored, and that the quenching works differently for different components.) 




Figure 2. Clear disagreement between mean-field model with cubic 
quenching formula (solid line) and the result from the helicity con- 
straint (dashed line). For comparison the helicity constraint formula 
with a 3t7 factor (instead of 2rf) is also shown as open symbols. 

One might think that a catastrophic quenching expression of the form 

1 — «b-B / B 2 q , as opposed to 1/(1 + obB /B 2 q ), for both a and r^o may also 
reproduce the resistively dominated growth of the field, but this is not the case. 
To show this we consider the steady state of the a 2 -dynamo in the one-mode 
approximation, so 



dlnB 2 



= (a k - r/tofe 2 )(l - a B B z /B 2 q ) - V k 2 , (4) 



dt 

which leads to the final field strength B^ n given by 



r2 /r2 i V k a o k ~ ^Tofc 2 _ . . 

a BB Rn t> =1 -r = —r K, 1, (5) 



where r/xo = r ?+ 7 ?to is the total (microscopic plus turbulent) magnetic diffusivity. 
Since we know from the simulations that B^ n /B eq = 0(1), Eq. (||) implies that 
as = 0(1) i which would then not reproduce the resistively dominated growth. 

In order to appreciate the difference to a quenching formula of the form 
1/(1 + a B B /B 2 q ), we write down the corresponding equation for that case: 

d* 1 + a B B 2 /B? 



eq 
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so 

= a k - r]tok 2 - r]k 2 (l + a B B 2 /Bl q ) = a k - r] T0 k 2 - r]k 2 a B B 2 / B 2 q , (7) 
which leads to a final field strength B^ n given by 

n R 2 /R 2 _ apk - i]Tok 2 

olbB^JB^ -p > 1, (8j 

so Bfi n /B eq = Oil) implies a B 1, as in Eq. (Q), which we know leads to 
resistively dominated growth. 

Although we cannot exclude that there may be other possible quenching 
formulae than Eq. (^), it is clear that it is easy to come up with other plausible 
choices which do not obey the helicity constraint ([l]). However, despite its suc- 
cess, (||) cannot be regarded as fundamental in that under other circumstances 
it may not describe the a-effect correctly. In the following we describe first an 
example where Eq. (|2|) gives a good prediction of what is then also verified in 
simulations, and then we turn to another example where Eq. (Q) does poorly. 

4. Open boundaries 

The main problem of why the saturation of the large scale dynamo progresses on 
a resistive timescale is, mathematically speaking, that there is no other term on 
the right hand side of the magnetic helicity equation. However, this changes if 
the volume under consideration is no longer periodic or closed, so that there can 
be a transport of magnetic helicity through the boundaries. However, what tends 
to happen is that the loss of helicity occurs mostly on large scales. Associated 
with this is a significant loss of large scale magnetic energy. Thus, although it is 
true that the time scale gets reduced, the saturation energy decreases, making it 
harder to reach significant field strengths of the large scale field. (The small scale 
field is unaffected by this.) In fact, looking at the field evolution on an absolute 
scale shows that the loss of magnetic helicity and energy simply terminates the 
growth at a lower level; see Fig. ||. 

One may speculate that lower saturation levels are still acceptable in many 
astrophysical settings, because if the dynamo is sufficiently supercritical (like in 
accretion discs and many stars, but possibly not in galaxies), the tendency to 
reduce the saturation level could be offset by a correspondingly stronger dynamo. 

5. The effects of shear 

The presence of shear is interesting because in that case oscillatory solutions can 
be expected. This is indeed what happens; see the butterfly diagram in Fig. |j. 
In this model the oscillation frequency is 27r/1000 ~ 0.006 (in nondimensional 
units). This value lies just between the ohmic and dynamic rates of change, 
r\k\ = 5 x 10~ 4 and 0.08, respectively. This suggests that the cycle period is not 
affected by the magnetic helicity constraint in the same way as the growth time 
of the dynamo, which is a resistive timescale. For a full account of this work see 
Brandenburg et al. (2000). 
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Figure 3. Evolution of magnetic energy (first panel: normalized to 
the final value, second panel: unsealed) for a mean-field model with 
catastrophic quenching and either open boundaries (solid lines) or pe- 
riodic boundaries (dashed lines). In the third panel we show a merid- 
ional cross-section of the toroidal magnetic field from a simulation with 
open boundaries in the z-direction at ±7r. Different shades represent 
different values of the toroidal field strength. The parameters of the 
mean-field model are ao = 2, r/ t o = 1, r] = 0.01, k = 1, so the kine- 
matic growth rate is A = a^k — rj^k 2 = 1, and with kf = 5 we have 
«B = \/{r]kkf) = 20. In the non-periodic case a vertical field boundary 
condition was adopted. 
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Figure 4. Butterfly (space-time) diagram of the toroidal mean field 
at two different x-positions where the sign of the shear is either negative 
or positive, giving rise to dynamo waves traveling in opposite direction. 
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6. Conclusions 

An important step has been reached in modelling dynamo action in astrophysics 
in that we are now beginning to bridge the gap between direct three-dimensional 
simulations on the one hand and the mean-field approach on the other. Many of 
the conclusions reached earlier in the framework of mean-field theory do apply, 
but there are also important restrictions. Most importantly, the alpha and eta 
coefficients have to obey some rather general restrictions resulting from helicity 
conservation. Further extensions to the theory, for example the inclusion of 
open boundaries, can straightforwardly be modelled. However, it seems that 
open boundaries do not enhance dynamo action, as one may have hoped, but 
instead they lead to a lower field strength at saturation, which is now however 
reached earlier on a dynamical timescale. In the presence of shear the helicity 
constraint still applies, but the modelling in terms of mean-field theory is not yet 
fully understood because anisotropics in the quenching become very important. 
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